Relaxation in driven integer quantum Hall edge states 
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A highly non-thermal electron distribution is generated when quantum Hall edge states originating from 
sources at different potentials meet at a quantum point contact. The relaxation of this distribution to a stationary 
form as a function of distance downstream from the contact has been observed in recent experiments [C. Al- 
timiras et al. Phys. Rev. Lett. 105, 056803 (2010)]. Here we present an exact treatment of a minimal model for 
the system at filling factor u=2, with results that account well for the observations. 
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Introduction. The importance of understanding non- 
equilibrium dynamics and relaxation in many-body quantum 
systems has been recognised since the early years of quan- 
tum mechanics Settings in which such problems are of 
high current interest include, among others, cold atomic gases 
^ and nanoscale electronic devices. jsl-lTill As a particular 
example, recent experiments |3] on quantum Hall (QH) edge 
states driven out of equilibrium at a quantum point contact 
(QPC) provide very detailed information on the approach to 
a steady state in an electron system that appears to be well- 
isolated from other degrees of freedom. In this paper we de- 
scribe the exact solution of a simple model for these experi- 
ments and compare our results with the measurements. 

In outline, the experiments we are concerned with |3] in- 
volve two sets of integer QH edge states, which meet at a 
QPC. When a bias voltage is applied to the QPC, tunneling 
between the edge states generates a non-equilibrium electron 
distribution. The form of this distribution in energy and its 
evolution as a function of distance downstream from the QPC 
are probed by monitoring the tunneling current from a point 
on the edge, through a quantum dot that has an isolated level 
of controllable energy. Close to the QPC, the measured distri- 
bution has two steps, reflecting the different energies of Fermi 
steps in each of the incident edges. With increasing distance 
from the QPC, the distribution relaxes to a single, broad step. 
The theoretical challenge presented by these observations is 
to understand and model this relaxation process. 

The relationship between these edge state experiments and 
other recent work on many-body quantum dynamics far from 
equilibrium has several aspects worth emphasising. First, in 
the context of QH edge states, these are the most recent of 
a series of striking observations of non-equilibrium effects in 
interferometers |i6J and in thermal transport |7J. They are also 
the equivalent for a ballistic system of earlier studies J^l of lo- 
cal distributions in diffusive wires. Second, the measurements 
stand apart from earlier theory |13| and experiment lfl4]l on 
non-equilibrium transport between fractional QH edge states, 
because they probe local distribution functions, rather than 
the global non-linear current-voltage characteristic. Third, 
and more broadly, the system studied is different in impor- 
tant ways from quantum impurity problems [o'-T?], as there is 
no impurity degree of freedom and interactions are not con- 



fined to an impurity site but instead operate in the ingoing 
and outgoing channels. Fourth, there is an analogy between 
edge state relaxation and cold atom experiments in the time 
domain |2], since distance from the QPC translates roughly as 
time, using the edge state velocity as a conversion factor In 
that sense the experiment we consider, probing relaxation as 
electrons propagate, is equivalent to a quantum quench jisll . 
in which time evolution is studied following a sudden change 
in the Hamiltonian. An important question in this context is 
whether a system thermalises at long times. Since integrabil- 
ity is an obstacle to thermalisation, it is noteworthy that the 
conventional model lll6ll of a QH edge state as a chiral Lut- 
tinger liquid is an integrable one. 

Any attempt to model theoretically the experiments of 
Ref. 3 starting from a chiral Luttinger liquid description faces 
an obvious difficulty, since interactions are most naturally de- 
scribed in terms of collective modes using bosonization, but 
tunnelling at the QPC is simple only in terms of fermionic 
variables. In pioneering work, two alternative approaches 
have been developed: one based on a Boltzmann-like equation 
for the electron distribution OlTIl : and the other using a phe- 
nomenological model for the plasmon distribution generated 
at the QPC |18]. From fis"], and from subsequent discussion 
of a quantum quench in an isolated QH edge (with an initial 
state chosen to emulate the effects of a QPC) by the present 
authors fl^, a physical picture has emerged, in which relax- 
ation is seen as a consequence of plasmon dispersion or the 
presence at = 2 of two plasmon modes with distinct veloc- 
ities. In summary, an electron that tunnels at the QPC can be 
viewed as a superposition of plasmons. Dispersion or multiple 
plasmon velocities cause such a wavepacket to broaden as it 
propagates. The lengthscale for relaxation of the electron dis- 
tribution downstream from the QPC is the distance at which 
the width of this wavepacket is comparable to the character- 
istic separation between tunnelling electrons. The argument 
identifies relevant scales but does not generate a prediction for 
the electron distribution and its dependence on distance from 
the QPC. One of our main aims here is to calculate this central 
quantity. 

An overview of our treatment is as follows: we make key 
use of the fact that, since experiments are at Landau level 
filling factor v=2, each QH edge carries two co-propagating 
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channels. Interactions mix excitations in different channels, 
and in the minimal model eigenmodes are charactised us- 
ing two velocities. Our approach combines refermionisation 
iEq. I2P1 with non-equilibrium bosonization methods I^ 22l - l26ll : 
by first bosonizing, then recombining bosons to form a new 
set of fermions, we transform the Hamiltonian for the inter- 
acting system into one for free particles. Under this trans- 
formation, the measured distribution (or more accurately, the 
spectral function probed by the tunnelling conductance) is ex- 
pressed as a free-fermion determinant. This determinant has 
a form similar to that appearing in the the theory [27, 28] of 
full counting statistics [FCS], and related quantities appear in 
a variety of other non-equilibrium problems 1 29] . Its numeri- 
cal evaluation can be done accurately and efficiently, yielding 
the results we present below. In addition, simple analytical ex- 
pressions can be found for some quantities, and we give one 
for the total energy in each channel as a function of distance 
from the QPC, thus providing a Hamiltonian derivation of a 
result obtained previously in the phenomenological treatment 
of Ref. 18. 

The work we set out here is related in a variety of ways to 
other studies of devices built from integer QH edge states. In 
particular, experiments that revealed striking non-equilibrium 
effects in Mach Zehnder interferometer s M have stimulated 



extensive theoretical research, E MM 

including calcu- 
lations ll33ll of interferometer dephasing based on the same 
bosonized model lildl for edge states at v = 2 with contact 
interactions that we adopt in the following. In the context 
of edge state relaxation a parallel development to the present 
paper, described in Ref. is based on an approximation (in- 
volving a factorisation of bosonic correlators) that is expected 
to be accurate sufficiently far from the QPC. Our solution re- 
covers such a factorisation, but only at distances large com- 
pared to the relaxation length for the model. 

Model. The experimental system is illustrated in the upper 
panels of Figs. [T] and |2] two alternative geometries arise, ac- 
cording to whether the tunnelling conductance is measured in 
the same or the opposite channel to that coupled by the QPC. 

We take a Hamiltonian with kinetic, interaction and tun- 
neling terms, in the form H ~ ffkin + ^^int + ^^tun- Using 
the labels r] — 1,2 to distinguish edges according to their 
source, and s =t, i to differentiate between the two channels 
on a given edge, the fermion creation operator at point x for 
channel 77, s is ^^^(x). It obeys the standard anticommutation 

relation {iipljg{x)ip,j's'{x')} = Sss'Sjjjj'S{x — x'). The density 

operator is p,-is{x) = ipjjg{x)i^jjs{x)- Taking all four channels 
7], s to have the same bare velocity v, and assuming a contact 
interaction of strength g between electrons in different chan- 
nels on the same edge, we have [ 16[] 

-ffkin = -«?iw^y" V'is(a;)3a;V'T,s(a;) da; (1) 
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and 



Hint = 2'Khg'^ J Pjj^{x)pjji{x) dx . (2) 




FIG. 1: (Color online) Top: sketch of the experimental geometry 
in which the tunnelling conductance is measured in the same chan- 
nel as that coupled by the QPC. Bottom: differential conductance 
calculated (lines) and measured |3] (symbols) in this geometry at in- 
dicated distances from QPC. Fits use T = 44 mK, V = 30.3 
(PeV = 8), and p — 0.545. Black lines: calculations (full) and data 
(dashed) for ^=0. 



Note that short-range intrachannel interactions can simply be 
absorbed into the value of v. Tunneling with amplitude ^qpc 
at the QPC between channels 1| and 2J, is described by 

i^tun = tQPci'l^ (0) V^ai (0) + h.c. . (3) 

A bias voltage V generates a chemical potential difference eV 
between incident electrons on edge 1 and those on edge 2. 

The observable of interest is the tunnelling conductance as 
a function of energy E and distance d > from the QPC, in 
channel 2, s. This is the Fourier transform of the correlator 
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-iHr/h 



(4) 



where the average is taken in the non-equilibrium steady state, 
with s^l in the geometry of Fig. [T|and s=t in that of Fig. |2] 
The tunnelling conductance is determined by a measurement 
of the current Iq£){E) through a quantum dot with a single 
level at energy E weakly coupled to the channel. With tun- 
nelling amplitude tD to the dot, this is lfl9ll 



(5) 



Results. We show in fss'] that the Hamiltonian H can be 
brought into a free particle form by introducing transformed 
fermion operators {x) with index a = A± , S± and eigen- 
mode velocities Va = v± = v ± g. Using these new operators 



[tQPC ^\40)^A-{0)+h.c.] 



(6) 
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Moreover the tunnelling conductance is obtained from 

Gs{d,T) oc (e--[±^^-K-)+-^^+('^.-)l) , 



(7) ^ 



where upper sign is taken in the first term of the exponent for 
s= t and the lower sign for s= J,. This expectation value is 
taken in the stationary scattering state of H specified by the 
temperature and chemical potentials of incident channels, and 
the operators 



d-\-v± r 



(8) 



count electrons within the spatial intervals d < y < d + v±t 
on the channels A±. 

Equations (|6]l - ([8]l form the central results of this paper: 
they express the observable of interest in an interacting, non- 
equilibrium system, in terms of the expectation value of a 
single-particle operator, evaluated as an average in a station- 
ary scattering state of a single-particle Hamiltonian. We now 
use these equations to discuss the physics of relaxation in this 
system and to calculate the tunnelling conductance. 

A simple physical picture of the relaxation process, and an 
identification of relevant length scales, follows from the form 
of Eq. (|7]i. This picture has the same content as our discussion 
(above and in Ref. of relaxation arising from two plasmon 
velocities, but is phrased in terms of the transformed fermions. 
Note first that the operators N± {d, r) count fermions that pass 
through the QPC in two time windows, both of duration r, but 
with a relative delay of d/vcs, where 



l/«cff = 1/v- - l/v+ = 2g/{v' + g^). 



(9) 



At d = the two windows exactly overlap. Then, for ex- 
ample, G|(0, r) acquires the same contribution from each 
fermion inside the window, regardless of tunneling. In con- 
sequence, and as expected, the electron distribution at c? = 
in the t channel is unaffected by tunneling, since the QPC 
couples I channels. (Conversely, it is clear that Gj. [d, r) is 
affected by tunnelling even for d 0, since it depends on 
the difference and not the sum of the fermion numbers in 
the two windows). Far downstream from the QPC, by con- 
trast, d/vcs ^ T, the two windows are widely separated in 
time, and the fermion numbers in each are uncorrelated. In 
this case, the right side of Eq. d?) is the product of indepen- 
dent factors, each with the form of a FCS generating function 
11271 12811 . At large d the tunnelling conductance is the same 
in both channels. It matches exactly the result obtained previ- 
ously by considering a quantum quench, with the impor- 
tant consequence that even the limiting distribution far from 
the QPC is non-thermal, although deviations from a Fermi- 
Dirac form are small, (l^ The velocity Vcs can be combined 
with voltage V or temperature T to define two lengths: 



ly = hvce/e\V\ and /t = hvcs/2TrkBT. 



(10) 




These set the scale for relaxation, and diverge in the non- 
interacting limit. 



FIG. 2: (Color online) Top: sketch of the experimental geometry in 
which the tunneling conductance is measured in the opposite channel 
to that coupled by the QPC. Bottom: differential conductance calcu- 
lated (lines) and measured 1 3] (symbols) in this geometry at indicated 
distances from the QPC. Fits use T = 44 mK, V = 45.5 fj,V, and 
p = 0.545. Black lines: calculations (full) and data (dashed) for 
V = 0. 



Detailed comparison of our results with the measurements 
of Ref. 3^ requires an evaluation of the tunnelling conductance 
at general d. For this we compute Gg {d, r) from Eq. (|7]i as 
summarised in BSll . The outcome is shown for the two al- 
ternative geometries in Figs. [T] and |2] The theoretical curves 
match the data well, except for the smallest value of d in 
Fig-H Four parameters enter the calculations: Vctf, T, V and 
the tunneling probability p at the QPC. In addition, one further 
parameter is needed to set the energy scale for the data: the 
gate voltage-to-energy lever arm rjQ (defined in Ref. 9, sup- 
plemental material). Of these parameters, only Vcff is uncon- 
strained by independent measurements. We arrive at the fits in 
Figs. [T]and|2]as follows. Wefix?7G = 0.062 and p~ 0.5 (con- 
sistent with the values determined in Refs.Qandj^) and obtain 
the bath temperature T = 44 mK (close to the experimental 
value of T = 40 mK) by fitting data at zero bias voltage to a 
thermal Fermi-Dirac distribution. We then compare numeri- 
cal results at large d with data measured in the | channel at the 
maximum distance from the QPC (30 /im). We obtain a good 
fit at /3eV ~ 8, and use this value for computations shown 
in Fig. [T] at all other d. Finally, we fix the interaction scale 
VcB by fitting data at one intermediate distance (4 /im). The 
theoretical predictions at other distances are then fully deter- 
mined. (Our specific value for p [0.545 rather than 0.500] is 
significant for the fit only at intermediate d and small energy.) 
Data for the geometry of Fig. |2] were taken using a different 
bias voltage from that in Fig.[T]and we calculate the theoretical 
curves shown in Fig. fusing the same values of j]q, T and Vcs 
as for Fig. [1] but scaling the value of (5eV by the ratio of ex- 
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perimental voltages. Both values of /SeV are 25% lower than 
the experimental ones, and this discrepancy is consistent with 
the 'missing energy' reported in Ref. Q Possible interpreta- 
tions of this discrepancy include loss of energy at the QPC 
or the existence of additional, neutral edge modes, illli3,[Iii 
or invoke 1 19, l^l the distinction in an interacting system be- 
tween the measured spectral function and the electron distri- 
bution in energy. The fact that we obtain good fits to all data 
except that in Fig.[T]for the smallest d suggests that, if energy 
is lost from the edge modes in the experimental system, this 
probably happens very close to the QPC. The most significant 
outcome from this fitting process is a value for the interaction 
parameter: Vctf — 6.5 x 10'*ms^^. A roughly comparable 
result (fcff ~ lO^ms^^) was obtained from the same data in 



Ref 1 1 8L but we believe our fitting procedure is more precise. 
Both determinations are in the range measured for edge state 
velocities in gated samples. [3] 

An additional quantity of interest is the energy density 
£"(6?) in channel s: it characterises with a single value the 
electron distribution |3]. Moreover, we can evaluate it ex- 
plicitly from Eq. ^ because (see Ref. Il9l) it is proportional 
to drG^{d,T)\r=o and G^{d,T) has an expansion for small 
T in cumulants of N±{d,T). Interactions cause energy ex- 
change between channels, so that £*((i) evolves with d. We 
obtain 



£^{d)=eT + — 



It 



/2 



Sin' 



d/2h 



sinh^ d/2h 



(11) 



Here ey = p{l — p)eV / AttHv is the excess energy in a single 
channel due to a double-step electron distribution generated 
by a QPC with tunneling probability p. Eq. (fTTl l. derived using 
a different approach in Ref. 18, shows that the energy densities 
interpolate between values just after the QPC of e^(0) — £t + 
£v and £^(0) = £t, and values at large d of £^(c?) = £t + 
£v/2. The relaxation length is set by the smaller of Zt and 
ly. This relaxation is oscillatory, but oscillations are strongly 
suppressed if eV ^ k-gT. 

Details of calculations are given in Supplementary Material 
1 38|] . We thank F. Pierre for extensive discussions and for pro- 
viding access to the data of Ref. (sl The work was supported 
in part by EPSRC grant EP/I032487/1. 
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SUPPLEMENTARY MATERIAL 



MODEL 

We take a Hamiltonian with kinetic, interaction and tun- 
nelling terms, in the form H = -ffkin + -f^int + -fftun- Us- 
ing the labels -q = 1,2 to distinguish edges according to 
their source, and s =tii to differentiate between the two 
channels on a given edge, the fermion creation operator at 
point X for channel rj^s is ^/'^^(x). The density operator is 

Pris [x) — V'j'js {x)4'ris {x)- Taking all four channels 77, s to have 
the same bare velocity v, and assuming a contact interaction 
of strength g between electrons in different channels on the 
same edge, we have llj] 



kin 



{x)dxipnsix) dx 



and 



Hi, 



2TThg 



(12) 



(13) 



Tunneling with amplitude igpc at the QPC between channels 
14, and 21 is described by 



Htun - iQPC^Ij,(0)V^2i(0) 



h.c. 



(14) 



A bias voltage V generates a chemical potential difference eV 
between incident electrons on edge 1 and those on edge 2. 

The observable of interest is the tunnelling conductance as 
a function of energy E and distance d > from the QPC, in 
channel 2, s. This is the Fourier transform of the correlator 



where, as usual, a is a short-distance cut-off, iV^s is the 
fermion number operator, and L is the length of the edge. 
After bosonization, the combination ifki„ + Hint is diag- 
onalized by a rotation to new bosonic fields Xa, where 

(XS+ XA_ XA+ Xs-V = U{^if $21 <i>2tV and 

/I 1 1 1^ 

1-1 1-1 
1 1-1-1 

VI -1-1 1/ 

Next we use the fields Xa to define new fermion operators^ 
\^ ^c(. Crucially, besides diagonaUsing iJkin+^int, the trans- 
formations ensure that Htun remains a single-particle operator 
when expressed in terms of ^a- Specifically, after bosoniza- 
tion 



(18) 



(19) 



while under rotation 0u(O) - 02^(0) = XA+{0) - Xa_(0). 
Since XA± (0) appear in this expression with unit coefficients, 
we can introduce new Klein factors Fn and new fermion fields 



Fr. 



to obtain the expression for H displayed 
in Eq. (6) of the main text. The required transformation of 
Klein factors has been given previously in the context of a 
two-channel Kondo model: |5] since fermion number opera- 
tors should transform following Eq. ( fTSl l. we require 



- fI^Fii, 



Fs_F 



t 

A- 



(20) 



Gs{d, t) = {e'"^l^i}l^{d)e''^"''/^i}2s{d)), (15) This implies that the combination appearing in Htun has the 



^2s 

where the average is taken in the non-equilibrium steady state. 
The tunnelling conductance is determined by a measurement 
of the current Iqu{E) through a quantum dot with a single 
level at energy E weakly coupled to the channel. With tun- 
nelUng amplitude to to the dot, this is 12] 



iQoiE) 



ft2 



Gs{d, T)e-'^^/^ dr. 



(16) 



REFERMIONIZATION 

The steps required to derive Eqns. (6) - (8) of the main text 
are as follows. First we bosonize the Hamiltonian in the stan- 
dard way [3], introducing bosonic fields (f>rjs{x) with com- 
mutation relations [(pjfs{x) , dyCflrjs' {x')] = —2TTi6jjjj'6ss'S{x — 

x'), and Klein factors Fjjs, and representing the fermion oper- 
ators via the relation 



ij^r,s{x) = (27ra) 



-1/2 F 



N„ 



-i(j>r,s (x) 



transformation F^^F2^ — 



f\^ Fa^ ■ Note that a unit change 
in the occupation number of one of the new fermions results 
in changes of one half for the occupation numbers of the orig- 
inal fermions. Physical states in the new basis must therefore 
satisfy certain selection rules, to ensure that fermion occupa- 
tion numbers in the original basis are integer These selection 
rules are set out in Ref. i5j. For our purposes the key point is 
that -fftun does not connect the physical and unphysical sec- 
tors, because the new fermion operators appear in it in pairs. 

We can also transform in this way the operators that are re- 
quired to generate an incident state with a density difference 

andF2;^2t = 



between channels. Since = F\^Fg 

Fa^ Fg^ , a bias voltage V between channels 1 and 2 is rep- 
resented by setting the chemical potential to be eV in channel 
and zero in channels and S±. 
The task now is to evaluate the correlation function 
Gs{d,t). For this it is convenient to work in the interac- 
tion representation, with -fftun as the 'interaction' and Hq = 
H]iin + Hint, using a superscript / to indicate operators in this 

Time evolution of 



(17) representation: A^{t) = e'^°^Ae~ 



iHnt 
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the new bosonic and fermionic fields fields is simple in this 
picture: we have 

XA±{x,t) = XA±{x-V±t), 

and similarly for Xg_^_ {x, t) and ^gj^ (x, t). 

The time evolution operator in the interaction representa- 
tion is 



=exp 



4'u„(T)dr 



(21) 



(the usual time-ordering is not required here because 
iJ/mj(t2)] = for all ti and ^2)- Scattering states 
are generated by the action of S^{t): Eq. (fTSl l takes the form 
Gs{d,T) = {Q)o, where 

Q^[5^(r)]t[v;L(d,r)]t5^(r)[5^(0)]t4^,(d, 0)5^(0) 

and (. . .)o denotes a conventional thermal average, with 
Hamiltonian Hq and chemical potentials iii and /i2 on the two 
edges, defined by 



(...)o = Z-iTr|e^''(-^«-^i 



Ni~p.2N2) 



}• 



(22) 



where Z = Trje-'^^^f^^'i^i-^^^*^)} and A^i = Ni^ + Ni^, 
the number operator for edge 1 (and correspondingly for N2)- 
We now describe how the quantity Q may be simplified. 
First, by a straightforward though lengthy calculation one can 
show that 

[HLn{ti),i'L{d,t2)] = for h i [t2-d/v+M-d/v4 . 

Next, we insert [5'^(oo)]^S'^(oo) between the factors of S^{t) 
and [5^(0)]^ in our expression for Q. We then commute 
S\t)[S\oo)\^ to the left, and 5^(oo)[S'^(0)]t to the right, 
to obtain 

g - [S\^)]'^[^,{d,T)]m^{d,Q)S'{^) . (23) 
Bosonizing using Eq. (flTl l we have 

Rotating to the new fields and using 

2-K - 
Xi (rf, 0) - xi (d, r) + — [Ni (r) - Ni (0)] = 27rAA„ {d, r) 

with Ma {d, t) defined in Eq. (8) of the main text, and setting 
/i2 = so that Fj^{T) — F2S, we find 



where the signs ± correspond to the choices s =t, Since 
(t) transforms only the channels A± and not S±, we can 
extract a normalisation factor to arrive at 



G,(d,r) =Go(r) 



(24) 



where (. . 

Go(r) 



; ([S'^(oo)]t ...S'^(oo))o and 

{i^lid,r)^2s{d,0))o 



1 



sinh2 



(— w_T + ia)] 



(25) 



A final step is to consider the effect of the evolution opera- 
tors [^^(cx))]^ and {00) on the basis states in which the ex- 
pectation value is calculated. As these states are generated by 
the action of the operators ^^^±(2;) on the vacuum, we must 
find how these operators transform. We do this by solving the 
Schrodinger equation with the Hamiltonian of Eq. (6) of the 
main text. The solution involves the scattering amplitudes at 
the QPC: introducing 6 = tqpc jh^Jv^AjZ, the tunnelhng and 
reflection probabilities are sin^ Q and cos^ 6*, while 

cos6'\E'A+(a;) - isin6' (26) 

and 

[^^(oo)]t*^_(a;)^^(oo) = 

cos 6* (a;) - j sin 61 *^^(w+a;/w_) . (27) 

Hence the average (...) in a scattering state is evaluated by 
combining the rotation between channels defined in Eqns. ( 
and ( |27] | with the thermal average (. . .)o. 



NUMERICAL EVALUATION OF THE TUNNELLING 
CONDUCTANCE 

In the following we outline the methods we use for numer- 
ical evaluation of the correlation function defined in Eq. (7) 
of the main text, and hence the differential tunnelling conduc- 
tance shown in Figs. 1 and 2 of the main text. The central 
problem is to evaluate the normalised expectation value ap- 
pearing on the right hand side of Eq. (l24l i. which has the form 



(28) 



with X = f.-iA^A^{.d,r)±HA-{d,r)\_ jjjg jjjj^jj of large dis- 
tance d from the QPC, contributions to the two time windows 
appearing in the exponent in X are independent and it is con- 
venient to use the approach described in Ref. At finite d 
there is no such factorisation and we resort instead to a method 
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similar to one described in appropriately modified to 

take account of finite temperature. 

(i) Long distance limit. At large d the expectation value 
factorizes into a product of two functions, each of which count 
number of particles in a fixed time window, so that 

with 

X±(5,r) = (e'^^±('^'-))„orm. (29) 

The functions x± t) have a FCS form and can be written 
in terms of determinants for non-interacting fermions. One 
hasil 

{e'SU^'^d^r)'^ = det{[l - F(e-'* - l)n±{e)P]} . (30) 

Here n± (e) is the corresponding electron energy distribution 
in a given channel at finite temperature and bias voltage after 
the action of 5^(oo) (i.e. a double-step) while P is a projec- 
tion operator that is diagonal in the time domain, having the 
action on a time-dependent function y{t): Py{t) = y{t) if 
t e [0, r] and Py{t) — otherwise. Using the regularization 
procedure proposed in Ref. we can write the determinant 
dSOl l in the form 

(e^^*±M)=det[/(t,-i,)], (31) 
where the function f{t) is defined as a Fourier transform of 

/(e) = [l-n±(e)(e-'*-l)]e-'ii, (32) 

which is periodic in the domain [—A, A], with A a high energy 
cutoff on single-particle states and tj = jV/A. Note that the 
phase factor e~*5f appearing here is crucial, since without it 
one would wrongly obtain a result periodic in S. 

The normalisation {e^^^^^'^''^^)Q is obtained from similar 
expressions in which the non-equilibrium distribution 
is replaced by a thermal one, and in practice we evaluate di- 
rectly the ratio x±i^i''')- We check numerical convergence 
by changing the value of the cutoff A. The size of the matri- 
ces required grows Unearly with r, the largest necessary being 
2000 X 2000. Some results obtained using this approach are 
shown in Figs. [3]and|4] The numerical calculations can be 
tested at large r by comparison with the asymptotic analytic 
results derived in Ref. 6. as shown in Fig. [3] the agreement 
is excellent. An incidental by-product of our calculations is 
the discovery of an interesting new feature of the asymptotic 
behaviour According to Ref. ^61 this is exponential in r, with 
a rate that diverges for 6 — ir and p = 1/2. In fact we find 
numerically that decay at these parameter values is not expo- 
nential but Gaussian in t. It would be interesting to search for 
an analytical derivation of that form. 

(ii) Finite distance. At finite distance d the correlation func- 
tion Gs {d, t) is not simply a product of independent contribu- 
tions from the two channels A±, and so cannot be expressed 




eVT/27T V 

FIG. 3: The quantity |x± (tt, t) p [see Eq. ( I29l >1 as a function of time 
T for different values of parameter j3eV at tunneling probability p — 
0.5 Thick solid line: jSeV = 2; dashed line: l3eV — 5; dot-dashed 
line: ^eV = 10; circles: fieV — 100. Thin solid line corresponds 
to asymptotic behaviour at large l3eV, which is represented within 
numerical errors by the function exp[— ^7r^r'^(fcBTds)^/fi^]. 




FIG. 4: (Color online). Comparison of the asymptotic behaviour of 
|x±(7i", t)!^ [see Eq. from Ref'^, (solid lines) with results ob- 
tained numerically using Eq. l l31t (dashed lines), at tunneling proba- 
bility p = 0.5 for different values of 5. From top right to bottom left: 
red lines S = n/2, green lines {S — Stt/I), blue lines (S = it — 0.1) 
and violet lines (5 = tt — 0.01). 

in terms of quantities familiar from the theory of FCS. We re- 
quire instead a different numerical approach. Using Levitov's 
determinant formula [7], we have 

= det[l-n(e)[[^(oo)]te-"[±^-^-('^'")+^^+('*'")l5^(oo)-l]. 

Here the determinant is in the two-channel Fock space and 
n{e) is a Fermi-Dirac distribution in each channel, but with 
two distinct chemical potentials. To evaluate this determinant 
we express AOii {d, r) as bilinears in the fermion creation 
and annihilation operators ^]^_^_{x) and ^a±{x). We eval- 
uate this determinant in a basis of eigenstates of if kin, con- 
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1 2 3 4 5 6 
eVT/2Tr V 

FIG. 5: Time dependence of the normalised correlator 

^g-i^KA+(d,r)-jV^_(d,r)]^^^^^^ j-j^j. jjjg channel as that 

coupled by the QPC. Curves are for tunnelling probability p — 1/2, 
eV/kBT — 5 and different distances d from the QPC: (thin solid 
line) d/lv — 0; (dashed line) d/lv — 2.5; (dotted line) d/lv = 5; 
(dot-dashed line) d/lv = 10; and (thick solid line) x/lv = 20. 

sidering edges of finite length with periodic boundary condi- 
tions and imposing energy cut-offs to obtain a matrix of finite 
size. For adequate convergence with increasing system size, 
we find that it is necessary to scale the lengths L± of the two 
channels according to the velocity, setting v+/L^ = 
With this choice, a basis of a few thousand states is sufficient 
to obtain the results we present here. 

We show representative results for the normalised correla- 
tor (e-*''[-^-^+(''^^)±^'*-('*^^)l)norm at finite d in Figs. Hand 
|6] The most distinctive feature is an oscillatory time depen- 
dence of the correlator in the channel that is coupled by the 
QPC, which at d = and T = can be obtained exactly as 




eW/27r ■() 



FIG. 6: Time dependence of the normalised correlator 

(^^-i^WA+(d,r)+AfA-(d,T)]^^^^^ fpj. jjjg channel not coupled 

by the QPC. Other parameters as for Fig. [5] 
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